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The differential cross section for radiative capture of protons by deuterons is calculated using 
different realistic NN interactions. We compare our results with the available experimental data 
below Ex = 20 MeV. Excellent agreement is found when taking into account meson exchange 
currents, dipole and quadrupole contributions, and the full initial state interaction. There is only a 
small difference between the magnitudes of the cross sections for the different potentials considered. 
The angular distributions, however, are practically potential independent. 


PACS numbers(s): 21.45.+v, 25.40.Lw, 25.20.-x, 27.10.4-h 


The radiative capture of protons by deuterons and the 
inverse reaction, the photodisintegration of ^He, have 
been investigated experimentally and theoretically over 
the last decades with quite some interest. Despite the 
various corresponding investigations, the theory is only 
in rough agreement with experiment, and there are incon¬ 
sistencies between the data up to 30% in the magnitudes 
of the cross sections. The experimental results by Belt et 
al. and King et al. ||,P are in good agreement. Those 
by Matthews et al. 0 and Skopik et al. Q agree in the 
angular distributions, but disagree in the magnitudes of 
the cross sections. This indicates a calibration problem 
of the measurements. 

From the theoretical side several attempts have been 
made to describe the cross sections in this energy region. 
In the early calculations by Barbour et al. ||] phenomeno¬ 
logical interactions were used. It was shown that the final 
state interaction is quite important, and that the E2 con¬ 
tributions in the electromagnetic interaction are needed 
in the differential cross section. In the calculations by 
Gibson and Lehman [Q a more realistic Yamaguchi in¬ 
teraction, but only the El components were employed. 
King et al. j|] performed an effective two-body, direct 
capture calculation with the initial state being treated as 
a plane wave, or as a scattering state generated from an 
optical potential. Fonseca and Lehman Q calculated the 
polarization observables Ayy and at the excitation 
energy = 14.75 MeV including only the El interac¬ 
tion. A calculation at Ex = 15 MeV based on realistic 
interactions and both, the El and E2 contributions has 
been done by Ishikawa and Sasakawa j|]. Another cal¬ 
culation of Ayy in this energy region is by Jourdan et 
al. [^. It was found in all these investigations that T 20 
is independent of the deuteron and the helium D-state 
probability, whereas Ayy shows a weak dependence on 
these quantities. 

Very-low-energy radiative capture processes are of con¬ 
siderable astrophysical relevance. The p-d radiative cap¬ 


ture, which at such energies is almost entirely a magnetic 
dipole (Ml) transition, was studied in plane wave (Born) 
approximation by Friar et al. [ pT| . In these investigations 
the authors employed their configuration-space Faddeev 
calculations of the helium wave function, with inclusion 
of three-body forces and pion exchange currents. Var¬ 
ious trends, e.g., the correlation between cross sections 
and helium binding energies, and their potential depen¬ 
dence were pointed out. More recently a rather detailed 
investigation of such processes has been performed by Vi- 
viani et al. |l^. Their calculations employed the quite 
accurate three-nucleon bound- and continuum states ob¬ 
tained in the variational pair-correlated hyperspherical 
method, developed, tested and applied over years by this 
group. 

In Refs. ||^,Q we have treated the ^He photodisinte¬ 
gration and the inverse radiative capture process within 
the integral equation approach discussed below. These 
calculations were based on the Paris, Bonn A, and Bonn 
B potentials in Ernst-Shakin-Thaler (EST) representa¬ 
tion: PEST, BAEST, BBEST lT|-0. We have demon¬ 
strated in particular the role of E2 contributions, meson 
exchange currents, and higher partial waves at E^ = 12 
MeV and E^ = 15 MeV. The sensitivity against the un¬ 
derlying potentials, moreover, was pointed out. In the 
present paper we extend these investigations and com¬ 
pare our calculations with all sufficiently accurate data 
below Ex = 20 MeV. 

The Alt-Grassberger-Sandhas (AGS) equations are 
well known to go over into effective two-body Lippmann- 
Schwinger equations |^| when representing the input 
two-body T-operators in separable form. The proton- 
deuteron scattering amplitude, thus, is determined by 

^(q,q") = v(q,q") -k J d^q' v(q,q') 5o(q') "^(q'. q")- 

( 1 ) 

Applying the same technique to the ^He photodisintegra- 
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FIG. 1. Total cross section for the capture of protons by 
deuterons. The data are from [1-5]. 


tion process, an integral equation of rather similar struc¬ 
ture is obtained for the corresponding amplitude Q|, 

■Ad(q) = 6(q)-h Jd^q'Goiq') Miq'). (2) 

In both equations the kernel is given by an effective 
proton-deuteron potential V and an effective free Green 
function ^o- However, in Eq. (|) the inhomogeneity 
of Eq. © is replaced by an off-shell extension of the 
^He photodisintegration amplitude in plane-wave (Born) 
approximation, 

B{q) = {q\{i^d\Hem\ipne)- (3) 

Here, [V'He) and \tpd) are the ^He and deuteron states, |q) 
is the relative momentum state of the proton. Hem de¬ 
notes the electromagnetic operator. In other words, with 
this replacement any working program for p-d scattering, 
based on separable representations or expansions of the 
two-body potential, can immediately be applied to cal¬ 
culating the full ^He photodisintegration amplitude with 
inclusion of the final-state interaction. 

The cross section for the p-d capture process is ob¬ 
tained from the corresponding photodisintegration ex¬ 
pression by using the principle of detailed balance 

_ 3 dcr“P 

dn “ 2 ^ dn ■ ^ ’ 

Here, k and Q are the momenta of the proton and 
the photon, respectively. In the present treatment no 
Coulomb forces have been taken into account. The ma¬ 
trix element for p-d capture differs from the corre¬ 
sponding n-d expression only in its isospin content. 

The results presented in this paper are obtained by 
employing the PEST, BAEST and BBEST potentials as 
input | p^ , p^ , however, with an improved parameteriza¬ 
tion by Haidenbauer . The high quality of this input 
has been demonstrated in bound-state and scattering cal¬ 
culations iIitI , ^ , ^ . 



FIG. 2. Angular distribution of the cross section for 
p-d capture at = 6.0 MeV. The data are from [2,3]. 


The electromagnetic operator relevant in the total 
cross section is, at the low energies considered, essen¬ 
tially a dipole operator. In the differential cross section 
we have to include also the quadrupole operator. Accord¬ 
ing to Siegert’s theorem Q, these operators are given by 




i=l 


and 


£f(2) 


E: 


2 3 


em \/20 ^ ^ G ^2x{'dij Pi') ^ 
^ i—1 


( 5 ) 

( 6 ) 


where E-^ denotes the photon energy, the nucleon co¬ 
ordinates, Ci the electric charges, and A = ±1 the polar¬ 
ization of the photon. 

Our method for determining the final state, i.e. the 
^He wave function, is described in Refs. 25|. In the 
calculation of the Faddeev components the total angular 
momentum j of the two-body potential was restricted to 
J <2, while in the full state all partial waves with j < 4 
(34 channels) have been taken into account. With this 
number of channels a converged calculation was achieved, 
incorporating 99.8% of the wave functions. 



FIG. 3. The angular distribution coefficient as func¬ 
tion of Ex- The data are from [1-5]. 
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FIG. 4. Differential cross section for p-d capture for energies from near threshold up to 16 MeV. The data are from [1-4]. 
The data set by Matthews et al. [4] has been renormalized with the Aq from King et al. [2,3]. 


Details concerning their high quality are given in [p^ . 
For the initial state all partial waves with j < 2 have 
been included in order to get a converged calculation of 
the cross section [|l3|Jl4| . 

Usually the differential cross section is expanded in 
terms of Legendre polynomials 

(t( 6») = Ao ^1-b (cos 0)^ . (7) 

The total cross section is obtained by integrating over the 
angle 9 between the incoming photon and the outgoing 
proton 

tr = dTr^Q. (8) 


Figure 1 shows the total cross sections for the Paris, the 
Bonn A and Bonn B potentials compared to the experi¬ 
mental data l-i- There is a small potential dependence 
of a and, hence, of Ag similar to the one observed in the 
corresponding photoprocess [p^pdf . In view of the error- 
bars, the experimental data by Belt et al. [Q and King et 
al. 1^,^ are reproduced for all potentials with the same 
quality. Those by Matthews et al. 0 and by Skopik et al. 
Q are not described by the theoretical curves. 

Figure 2 shows the angular distribution of the differen¬ 
tial cross section, i.e., the ratio of a{d) and the coefficient 
Aq compared to the experimental data ^ . This distribu¬ 
tion is evidently potential independent. In other words, 
its shape shows no correlation with the helium binding 
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energy, or the D-state probability of the ^He wave func¬ 
tion. 

Figure 3 shows the angular distribution coefficient 02 of 
the expansion (|^) compared to the coefficients extracted 
from experiment In accordance with Figure 2 there 

is almost no potential dependence, i.e. no dependence on 
the three-body binding energy and the H-state probabil¬ 
ity, although this probability varies for the three poten¬ 
tials considered between 6 to 8% [ p5[ . 

Figure 4 shows the differential cross sections obtained 
for these potentials at various energies compared to the 
experimental data. Due to the slight potential depen¬ 
dence of the total cross section and, thus, of Aq, the 
magnitudes of the curves differ correspondingly. In all 
cases there is good agreement between theory and exper¬ 
iment. As pointed out in this agreement can only 

be achieved by taking onto account El and E2 contribu¬ 
tions of the electromagnetic interaction, meson exchange 
currents, and higher partial waves in the potential and 
in the three-body wave function. It should be mentioned 
that for increasing energies the peak is slightly shifted 
to the right-hand side, because of a smaller El and a 
somewhat higher E2 contribution. Note that, due to the 
missing E1-E2 interference term, the quadrupole contri¬ 
bution is practically negligible in the total cross section. 

In [|^,Q we have shown that for different potentials 
the low-energy peak heights of the ^He photodisintegra¬ 
tion cross sections are strictly correlated with the cor¬ 
responding ^He binding energies, and with the number 
of partial waves included. The magnitude of the present 
radiative capture process, i.e., the constant Aq, appears 
to be similarly fixed by the three-body binding energy. 
In other words, at the energies discussed, the radiative 
capture cross section does not represent an additional 
observable for testing different potentials. 
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